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ABSTRACT 

A total of six deep exposures (using AOT CAM01 with a 6” PFOV) through the 
ISOCAM LW10 filter (IRAS Band 1, i.e., 12 pm) were obtained on a ~15 arcminute 
square field centered on the ecliptic plane. Point sources were extracted using the 
technique described by Desert, et al. (1999, A&A 342, 363). Two known asteroids 
appear in these frames and 20 sources moving with velocities appropriate for main belt 
asteroids are present. Most of the asteroids detected have flux densities less than 1 
mJy, i.e., between 150 and 350 times fainter than any of the asteroids observed by 
IRAS (Tedesco . et al., 2002a, Astron J., submitted). These data provide the first direct 
measurement of the 1 2 pm sky-plane density for asteroids on the ecliptic equator. 

The median zodiacal foreground, as measured by ISOCAM during this survey, is 
found to be 22.1 ±1 .5 mJy per pixel, i.e., 26.2 ± 1 .7 MJy/sr. 

The results presented here imply that the actual number of kilometer-sized asteroids 
is significantly greater than previously believed and in reasonable agreement with the 
Statistical Asteroid Model (Tedesco . et al., 2002b, Astron J., to be submitted.). 

1. INTRODUCTION 

Most main belt asteroids are found between 2.2 and 3.4 AU from the Sun and at 
ecliptic latitudes less than 20 degrees. Except for the largest asteroids, the actual 
number above a given size is poorly known. For example, estimates of the number of 

main-belt asteroids with diameters larger than 1 km range from 3x1 0 4 to IxlO 7 
(Farniella and Davis, 1994). 

The asteroid size distribution is important because it provides constraints on models 
of the original size distribution of the planetesimals formed in the inner solar system and 
their subsequent evolution. It is also an important datum in modeling the numerical size 
of the population of near-Earth asteroids and accounting for their evolution from the 
main belt into Earth-orbit-crossing orbits. 


’Based on observations with the Infrared Space Observatory (ISO), an ESA project with instruments 
funded by ESA Member States (especially the PI countries: France, Germany, the Netherlands, and the 
United Kingdom) with the participation of ISAS and NASA. 
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A typical 1 km main belt asteroid has a V magnitude between 21 and 23. It would be 
a straightforward program to survey, at visible wavelengths, all asteroids in given 
regions of the sky brighter than this limit. However, because, for any given distance, 
visual surveys are biased in favor of discovering higher-albedo asteroids, magnitude 
data alone cannot be used to accurately derive asteroid diameters . This is because the 
absolute brightness of an asteroid depends upon its cross section and albedo, and 
asteroid albedos span a range of at least a factor of 25 (0.02 to >0.5). Moreover, there 
may be systematic trends of albedo with size (Tedesco . 1994). 

Observing thermal emission permits us to obtain an accurate distribution of asteroid 
diameters because, unlike the linear dependence with albedo at visual wavelengths, the 
infrared flux is only weakly dependant on the geometric albedo. For example, on 26 
March 2001 the 100 km main belt asteroid 50 Virginia was at a solar elongation of 110°, 
a typical elongation for space-based infrared observations. Virginia’s visual magnitude 
at this time, given its SIMPS (Tedesco. et al., 2002a) diameter of 99.82 km, would be 
14.9 if its visual geometric albedo were 0.02 and 18.4 if its visual geometric albedo 
were 0.50, a difference of 3.5 mag. The 12.0 pm magnitudes under these same 
conditions, would be 2.16 and 2.46, respectively, or a difference of only 0.3 mag. 
Furthermore, the lower albedo actually results in a slightly higher 12.0 pm brightness 
because in this case the asteroid’s temperature would be higher. Thus, an IR survey is 
slightly biased in favor of discovering lower-albedo asteroids. 

To date there have been three space-based infrared surveys in which asteroids 
have been incidentally observed: The Infrared Astronomical Satellite (IRAS) , the 
Midcourse Space Experiment (MSX, Mill et al., 1994), and the Infrared Space 
Observatory (ISO) spacecraft, reported on here. For a description of the ISO mission 
see Kessler et al. (1996) and for details on the ISOCAM instrument see Cesarskv et al. 
(1996). 

Results on IRAS asteroids are given in Tedesco et al. (1992 . 2002a) and those on 
asteroids observed by MSX by Tedesco et al. (2002c). 

The IRAS asteroid survey is severely incomplete at low flux levels, i.e., below about 
1 Jy, because IRAS could only detect an asteroid in its survey mode if a known orbit 
was available. Thus, although IRAS observed at infrared wavelengths, it was limited by 
the (albedo-biased) visual surveys in which asteroids are discovered and from which 
data their orbits are calculated. (IRAS knowingly discovered no main-belt asteroids due 
primarily to the poor spatial resolution of its detectors.) 

IRAS and MSX serendipitously observed numerous asteroids in the course of their 
nominal missions. However, due to the way in which their observations were conducted, 
only asteroids with known orbits were identified with the infrared sources these 
spacecraft detected. IRAS observed -95% of the sky and MSX about 10%. Although 
the faintest asteroids detected in these surveys have flux densities of about 150 mJy, 
they are in no way complete to this flux level. The ISO asteroid “survey” (discussed 
below) observed about 0.125 sq deg of sky to a completeness limit of -0.6 mJy. 

2. THE ISO DEEP ASTEROID SURVEY - IDAS 

The goal of this survey was to cover the maximum area of sky to the faintest flux 
limit possible under the constraints imposed by the zodiacal background and the 
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available observing time. The field was selected to be in the ecliptic plane, near the 
upper limit of the ISO solar elongation constraint (i.e., near 106°), and located west of 
the Sun (to facilitate ground-based follow-up). In addition, the field was chosen to lie far 
from the Galactic plane and to contain no known IRAS sources or bright stars. The 
sensor used was ISO’s AOT CAM01 with a 6-arcsecond PFOV and using the ISOCAM 
LW10 filter (IRAS Band 1, i.e., 12 pm). 

Asteroids move and their flux may vary appreciably on time scales as short as 
minutes. Consequently, the exposure time was chosen to freeze asteroid motion on 
each sub-map, where a sub-map is a three-by-three arcminute area (the size of the 
ISOCAM array) in which each point was observed three times. Each sub-map consisted 
of a 30 second exposure sequence 2 at a fixed position followed by a step of one 
arcminute in ecliptic longitude or latitude where the exposure sequence was repeated 
and ending with another one arcminute step in ecliptic longitude or latitude where the 
exposure sequence was again repeated. The total time spent doing each point in a sub- 
map was 90 seconds. Because the apparent rate of motion for main belt asteroids 
(MBAs), under the observing geometry described above, is between 0 and 60 
arcseconds per hour, the maximum angular distance moved during the time required to 
obtain a sub-map is less than 1 ,5-arcsecond. However, each sub-map was sampled 
three times to create the complete map, and the times between successive sub-maps 
varied from 30 sec to 870 sec. Thus, the maximum distance a main belt asteroid would 
move between samples of a given point in the map is 14.5-arcseconds. 

Figure 1 is a schematic diagram of the map coverage. Each box is one arcminute on 
a side with North up and East to the right. The raster began with the 3’x3’ array located 
in the NE corner of the map, as indicated by the heavy lines around the nine cells in the 
upper right of the figure. One exposure sequence was made at this position and then 
the array was moved one arcminute (cell) west. Seventeen exposure sequences were 
made along a line of constant ecliptic latitude. This brought the array to the end of the 
first row. At this point, it stepped South one arcminute, made an exposure sequence at 
this position, and then made 16 one arcminute steps to the East to complete the 
second row. This process was then repeated until the center of the array had scanned 
17 rows. 

See the movie for an animated version of this figure speeded up by a factor of about 
60. As can be seen from the movie, or the numbers in the figure, each cell around the 
outer one arcminute of the map received no more than three exposure sequences and 
those in the one-arcminute border interior to this region no more than 6. All other cells 
in the map received 9 exposure sequences for a total of 180 seconds each. We refer to 
the region with 9 exposure sequences as the region of complete coverage. 

The 15x15 raster map was obtained in the same way but using 15, instead of 17 
steps. The complete coverage area is 1 5 sq arcminutes for the maps obtained with the 
17x17 raster and 13 sq arcminutes for those obtained with the 15x15 raster. 

Two maps as described above were made in June 1996 and another four in June 
1997. A total of 13.64 hours was expended in obtaining the observations presented 


2 A “clean" was performed at the start of each map. This flashed the array to remove the memory of the 
previous observation and required about 240 seconds. Next 25 stabilization frames were taken. The 
actual observations consisted of four 5-sec exposures at each point in the raster, plus an average of ten 
seconds to step to the next raster position. 
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herein. Maps 1, 2, 5, and 6 required 2.44 hours per map, while maps 3 and 4 had 
available 1.94 hours each. A consequence of the decision to keep the exposure per 
map point constant over all six maps was that maps 3 and 4 cover less area. 

The intention was to have maps 3 and 4 (the time for which was granted under a 
Supplemental Observing Request) made at least 36 hours after the end of the previous 
map pair. However, they were scheduled less than 12 hours after completion of the 
previous map and by the time the observing schedule was issued it was too late to 
reschedule them. 

The sample consists of six data sets that are now in the public domain 3 labeled as 
Target Dedicated Time (TDT) numbers: 21103003, 21103004, 57200101, 57200102, 
57200407, and 57200408 (corresponding, respectively, to map numbers 1,2,3, 4, 5, 
and 6 in Table 1). Thus, they were taken in pairs during two 24-hour ISO orbits (the first 
3 digits in the TDT) separated by approximately one year, on 15 June 1996 and 10 
June 1997 (i.e., on Julian days 2450249 and 2450609, respectively). 

Figure 2 shows the six images obtained after processing using the technique of 
Desert, et al. (1999), which is further described in Sec. 3. Figure 3 shows all point 
sources with SNR > 3.0 extracted from the ISOCAM maps shown in Figure 2. Squares 
outline the areas sampled nine times. The point size is proportional to the flux density, 
which ranges from 0.3 to 12.2 mJy 

3. ISOCAM DATA REDUCTION 

3.1 Observation characteristics 

A typical dataset consists of 1800 readouts, each with 5.1 seconds of integration, 
through the ISOCAM LW10 filter centered at 12 pm with a bandpass very similar to the 
IRAS 12 pm band. The lens wheel was on the LGe6 position providing a ratio of 6 
arcseconds per detector pixel (which is also close to the FWHM of the Airy pattern of 
ISO). The camera detector consists of 32 by 32 pixels, with one column (number 24, 
disconnected before launch) missing, providing a 3.2 by 3.2 arcmin instantaneous field 
of view. The total survey area was covered by making a raster with ISO at positions on 
a 17 by 17 grid with 60 arcsecond (10 pixel) steps and 60 arcsecond line separation. 
Each position was observed for four readouts, i.e., 20 seconds of integration time. With 

the survey redundancy (a factor 9), the total integration time per sky pixel is about 3 
minutes. The median zodiacal foreground, as measured by ISOCAM during this survey, 
is found to be 22.1 ± 1.5 mJy per pixel, i.e., 26.2 ± 1.7 MJy/sr (the error bar being the 
dispersion among the six surveys of the same area). 

3.2 Summary of data reduction 

The raw data consists of a cube (CISP files) of detector readouts (one every 36 
CAM time units i.e., 5.1 seconds) and an ISO pointing history (IIPH) file. The detailed 
data reduction procedure is described bv Desert et al. (1999). Here we give a summary 
along with the specific parameters that were used for the present datasets. First cosmic 


3 http://www.iso.vilspa.esa.es/ 
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rays are removed by a time line analysis of each pixel. Long duration glitches are also 
removed and a transient correction applied, using the method described by Coulais and 
Abergel (2000). 

The data timeline is then analyzed with a “triple-beam” linear algorithm that 
basically finds, for each camera pixel, the difference between the signal at one raster 
position and the average of the 2 adjacent position signals. The dispersion of this 
difference for different raster positions indicate the true pixel noise of the 
measurements, because, most of the time, it is uncontaminated by sources. Badly 
behaved pixel values (due to glitches and bad triple-beam X 2 ) are discarded by an 
adapted sigma clipping. We project the difference and dispersion on a final sky map, 
using neighbor pixel approximation, with a 2-arcsecond pixel size. A redundancy 
number is also obtained this way. The projection is done by coadding with an optimal 
weighing and a first order array distortion correction. We used the associations with the 
USNO optical catalog (version A2.0 - Monet et al., 1998) 4 5 to deduce the offset positions 
(up to 7-arcseconds in both directions) to apply to each dataset map (because of the 
so-called lens filter wheel jitter). 

The final map is then searched for point sources in a selected area where the 
redundancy is two or more. As explained by Desert et al. (1999), we iterate an 
algorithm where a candidate source (found with a top hat wavelet) is fit with a 9- 
arcsecond FWHM two-dimensional Gaussian (for the position and intensity) and the fit 
is removed. This algorithm allows measuring source fluxes near undefined pixels 
without underestimating the flux (as aperture photometry would do). The noise in the 
flux measurement is deduced from the noise map and the Gaussian least-square fitting 
algorithm. The absolute fluxes were deduced using the nominal ISOCAM internal unit to 
mJy conversion factor (/.e., by assuming that the factor has not changed with respect to 
pre-flight expectancy), and by applying a correction factor (1.52) to go from our fitted 
Gaussian beam flux to total point-spread-function integrated flux. 

In Table 2 s , we give the complete catalog of (527) sources that were detected at the 
£ 3 a level in any of the six maps. Column (1) is an identification number; columns (2) 
and (3) the J2000 RA and Dec; column (4) the flux density in band LW10, (5) the one- 
sigma uncertainty in the flux density; (6) the signal-to-noise ratio; (7), (8), and (9) are 
quality flags; (10) the Julian Date of the observation (an average of the, up to, nine 
measurements on each point source that are available); (11) a confusion flag; and (12) 
a code indicating whether the source was in the multiply-sampled region. Columns (13) 
through (19) provide data on sources found within 6-arcseconds of a USNO-A2.0 
catalog (Monet et al., 1998) visible source (as reduced by CDS-VizieR, http://cdsweb.u- 
strasbg.fr/CDS.html) . Column (13) gives the red magnitude from the USNO-A2.0 
catalog; column (14) the name from the USNO-A2.0 catalog; (15) the number of USNO- 
A2.0 sources associated with the ISO source; (16) the number from the USNO-A2.0 
catalog; (17) the distance from the USNO-A2.0 catalog source (95% are within 4- 
arcseconds); and (18) and (19) the distances in RA and Dec, respectively, from the 
USNO-A2.0 catalog source, rounded to the nearest arcsecond. 


4 As reduced by CDS-VizieR, http://cdsweb.u-strasbg.fr/CDS.html 

5 Table 2 is presented in its entirety in the electric edition of the Astronomical Journal. 
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A catalog of 63 inertial sources is given in Table 3. These are sources from Table 2 
that are seen at the same position in at least two maps. An average flux, ranging from 
0.34 to 10.8 mJy, and error is given along with the USNO-A2.0 catalog association. 
Some well-detected sources have no optical counterparts. These are probably external 
galaxies or very slow-moving asteroids. 

The complete (i.e., multiply observed) survey area is 225 arcmin' 2 (0.0625 deg' 2 ). 
The densities of stars and galaxies (assuming all non-USNO-A2.0 sources are 
galaxies), respectively, are found to be 0.072 ± 0.017 and 0.045 ± 0.013 arcmin' 2 , for 12 
pm flux densities greater than 0.6 mJy (a value close to the 4a level where a is the 
median flux error in the complete area for one dataset). 

4. ASTEROID IDENTIFICATION 
4.1 IDAS asteroids 

Moving sources were searched for in those areas observed in common within a day 
of each other. Twenty objects were found to move significantly (i.e., by more than 6- 
arcseconds over a period £ two hours). 

Because asteroids are moving sources, the technique described in Sec. 3.2 for 
obtaining the flux from coadded images underestimates their flux. Thus, we derived a 
rate-of-motion dependent correction factor (FDCor) by offsetting the inertial sources by 
various amounts to simulate their motion and then performing the photometry as 
normally on the coadded map. This resulted in smaller flux values as a function of the 
amount offset, to which we fit a second order polynomial (shown in Figure 7), viz., 

FDCor = 1 .001 - 0.0034*x + 0.01 60*x A 2, (1 ) 

where x = RT*30*Rate; and RT = 17 (for a 17 by 17 raster) or 15 (for a 15 by 15 raster), 
30 is the time per sample, and Rate is the apparent rate of motion in arcsec per revisit 
interval. 

Table 2 contains the uncorrected flux densities and Table 4, which presents the data 
on the 20 sources identified as being asteroids on the basis of forming tracks with two 
or more sightings, gives FDCor and the corrected flux densities, ranging from 0.43 to 
5.7 mJy, for each sighting. 

All of the identified asteroids have SNR > 4. The 1996 field contains four tracks in 
which at least one sighting has a flux density (FD) > 1 mJy, and the 1997 field contains 
six such tracks. Each field contains ten tracks, at least one of which has FD > 0.6 mJy, 
the 4-o completeness limit. 

Combining the results from the two fields gives 5 ± 1 probable asteroids with FD > 1 
mJy and 10 + 1 with FD > 0.6 mJy. 

Normalizing these results, gives 80 + 16 asteroids with FD > 1 mJy per sq deg at the 
ecliptic plane, i.e., with diameters greater than about 1.4 km at mid-belt, and 160 ±16 
with FD > 0.6 mJy per sq deg (diameters greater than about 1.1 km at mid-belt). 

Singletons (i.e., a source detected only once) with flux densities above 0.6 mJy may 
also be present in these fields but they cannot be unambiguously identified using these 
data alone. 


8 


4.2 Known asteroids associated with IDAS sources 


We associated two asteroids found in the 1996 field with known asteroids. Details 
on these sources are given in Table 5, where column (1) gives, as two rows per 
observation, the IDAS asteroid number assigned in Table 4 in the first row and the 
associated asteroid’s designation in the second, column (2) gives the IDAS source 
number from Table 2, columns (3) and (4) give, respectively, the observed RA and Dec 
in the first row and the predicted RA and Dec, obtained using the Horizons software 
(Gioraini et al., 1996) 6 , in the second row, columns (5) and (6) the observed corrected 
flux density and SNR, respectively, and column (7) the UTC of the observation. 

Of the ten probable asteroids identified in this field, those associated with 1999 AQ23 
and 17971 1999 JZ50 are the brightest. The predicted V-band magnitudes of these two 
asteroids at the time of the ISO observations were 18.2 and 18.1, respectively. This 
means that 80% of the asteroids in the 1996 field have V > 18 and, for those with low 
albedos (0.02), the maximum V is about 25. 

Figure 6 shows the observed ISOCentric positions for sources 3 and 1008 (1999 
AQ23) and sources 4 and 1003 (17971 1999 JZ50) from the 1996 ISOCAM map, 
together with the ISOCentric positions for the two known asteroids. According to the 
Horizons documentation: “The database is updated almost daily with new objects and 
orbit solutions. Comet and asteroid orbits are integrated from initial conditions stored in 
the JPL-maintained DASTCOM database 7 .”. Due to the ephemeral nature of these 
orbital elements, we present, in Table 6, those used in the analysis described here. 

The ISOCAM coordinates, for the mean time of the extracted sources, are plotted in 
Figure 6 as triangles. Because the asteroid position is the mean from detections 
obtained over the ~18 minutes required to map an inertial point on the sky, the 
predicted positions are shown as a series of 19 positions at one-minute intervals 
centered on the mid-time of the local map. The small squares centered on the predicted 
position trails are 6” on a side, the size of an ISOCAM pixel used in this experiment, 
while the figure is -3’ on a side, the size of the ISOCAM array. 

1999 AQ23 moved a distance equivalent to the size of a pixel during the time 
required to complete the raster scan of its position, while 17971 1999 JZ50 moved 
about twice this distance. For the numbered asteroid, 17971 1999 JZ50, the difference 
between the observed ISOCAM positions and the predicted ephemeris positions are 
2.8” and 0.6”. For the unnumbered asteroid 1999 AQ23, the observed positions lead 8 
the predicted positions by about 11” (or 34 minutes in time). 

As noted in Sec 3.2, the astrometric accuracy of the IDAS positions is better than 4” 
for 95% of the sources associating with USNO-A2.0 sources. The accuracy for moving 
sources is undoubtedly less, but probably not by a factor of two to three. The formal 
accuracy of the asteroids’ predicted positions is less than 1” (based upon output from 


6 The "HORIZONS (http://ssd.jpl.nasa.gov/horizons.html) On-Line Ephemeris System" was created and is 
maintained by the Solar System Dynamics Group, Jet Propulsion Laboratory. 

7 ftp://ssd.jpl.nasa.gov/pub/ssd/Horizons_doc.ps - Version 2.80 June 14, 2000. 

8 The lead-time is the difference between the time the known asteroid is closest to the position of the ISO 
source minus the time of the ISO observation of that source. The track of 1999 AQ23, in ISOCentric 
coordinates, passes less than 1” from the observed positions. 
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the Lowell Observatory Asteroid Ephemeris (ASTEPH) Version 1 .5 at 
httD://asteroid.lowell.edu/cai-bin/koehn/asteph run on 18 September 2001). 

In spite of the relatively poor agreement in position for 1999 AQ23 we nevertheless 
present the albedos and diameters for both 1999 AQ23 and 17971 1999 JZ50 under 
the assumption that they are associated with the ISO sources indicated. 

The results are presented in Table 7, where column (1) identifies the asteroid; 
columns (2) through (5) give the absolute visual magnitude, H, heliocentric distance, H 
Dist, geocentric distance, G Dist, and solar phase angle, Phase, from the Horizons 
ephemerides for the mid-time of the ISO observations; column (6) the corrected 
observed mean ISO LW10 band (-IRAS) 12 pm flux density and its uncertainty; 
columns (7) and (8) the computed diameter and geometric albedo using the Standard 
Thermal Model (STM - Lebofskv et al., 1986), D-STM and Ph-STM, respectively; and 
columns (9) and (10) the computed diameter and geometric albedo using the Near- 
Earth Asteroid Thermal Model (NEATM - Harris . 1998), D-NEATM and p H -NEATM, 
respectively. 

The results from the photometry are ambiguous. Using the given values for H and 
the infrared fluxes, the derived geometric albedos, from either the STM or NEATM 
thermal model, are either unphysical (-0.9) or implausible (-0.65 to 0.71). 

Physically plausible albedos can be obtained, e.g., by assuming that H for each of 
these asteroids is 1 .0 mag higher than that published, or that H is 0.5 magnitude larger 
and the 12 pm flux density is 50% higher than reported here, etc. (see Table 7). 
Changes of this magnitude for H are not uncommon (see Tedesco . 2002a). And a 50% 
underestimate of the ISO flux density is also quite possible. Furthermore, both of these 
asteroids are located in the inner part of the asteroid belt where, at least for asteroids 
with diameters greater than -60 km, low albedo asteroids make up less than 10% of the 
population (Gradie and Tedesco, 1982). However, other than ruling out low albedos for 
these asteroids (because for these values of H and an albedo of 0.02 the predicted 
infrared flux densities range from 200 to 440 mJy, two orders of magnitude higher than 
observed) an accurate albedo cannot be determined. 

5. SUMMARY 

There are about 1 60 asteroids per sq deg at the ecliptic plane above the ISOCAM 
LW10 band (/.&, IRAS Band 1, 12 pm) detection threshold of about 0.6 mJy. This 
corresponds to diameters greater than about 1.1 km at mid-belt. For the fields observed 
in this experiment the faintest asteroid source extracted has a flux density of 0.432 ± 
0.085 mJy and SNR = 4.3. 

To put these results in perspective, note that IRAS’ limiting sensitivity was about 150 
mJy at 1 2 pm whereas most of the asteroids detected by ISO are between 1 50 and 350 
times fainter. 

The Statistical Asteroid Model (Tedesco et al. 2002b) was run twice on a four sq. 
deg field centered on the ISO field, once for the epoch of osculation of the model’s 
orbital elements (14 Oct 1998) and a second time (yielding the results in parentheses 
below) for the date of the June 1 997 ISO observations. 

A total of 1 ,063 (1 ,638), of the approximately two million asteroids in the Statistical 
Asteroid Model, were present in this four sq. deg field. Of these, 673 (852) had 
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predicted IRAS 12 pm flux densities greater than 0.6 mJy. Thus, the model gives 180 ± 
20 asteroids per square degree with 12 pm flux densities greater than 0.6 mJy, in 
reasonable agreement with the ISO results reported here. 

The Statistical Asteroid Model result is actually a lower limit because the model does 
not yet include the NEAs or asteroids beyond the Hilda group. More importantly, 
however, is the fact that it terminates abruptly at a diameter of 1 km. If smaller asteroids 
were included some of these would have 12 pm flux densities greater than 0.6 mJy if 
they were close to the Earth. Nevertheless, the ISO data imply that the actual number 
of kilometer-sized asteroids is in reasonable agreement with the Statistical Asteroid 
Model. 

SIRTF’s imagers will have ~5’x5’ FOV arrays with 1.2” pixels and sensitivities of 
-0.015 mJy at 8 pm and -0.37 mJy at 24 pm, about an order of magnitude more 
sensitive than the IDAS limit. The results presented here demonstrate that significant 
numbers of asteroids will be present in virtually all deep exposures taken near the 
ecliptic plane. Using the extreme size-frequency distributions and the IDAS asteroid 
sky-plane density implies that there will be between 3 and 30 asteroids in each limiting 
sensitivity SIRTF image. See Tedesco et al. (2002b) for additional discussion on this 
issue. 

Unfortunately, in spite of the high quality of the ISO data, little can be said regarding 
the diameters (and nothing regarding the albedos) of the unknown asteroids detected. 
This is because their orbits, and hence distances and phase angles, are unknown and 
cannot be reliably computed from the ISO positions alone. Reliable diameter 
determinations for these ISO asteroids will have to await their discovery. And, with V 
mags probably fainter than 22 for most, this is unlikely to happen in the near future. 

In order to fully exploit such space-based infrared data, orbital elements of the 
asteroids must be known. If albedos are to be obtained, then visual wavelength 
observations are required as well. The minimum requirements for obtaining asteroid 
diameters in the absence of supporting ground-based observations is that the space- 
based data must sample the asteroid’s thermal spectrum at a minimum of three 
wavelengths bracketing the peak emission and be taken at appropriate intervals and 
with astrometric accuracies sufficient to allow computation of an approximate orbit. 
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Tables 


Table 1. Map Field Centers 


Map 

RA J2000 
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Raster 
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Longitude 

Ecliptic 

Latitude 
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Table 2. All Point Sources Extracted From the Six IDAS Maps 


This table is available only as a machine-readable table in the on-line electronic edition. 
































Table 3. IDAS Field Inertial Point Sources 


ID RA 

0 339.502145” 


Dec 

.506749 


FD o SNR 


USNO Name 


Associations 


339, 

.535786 

339, 

.540586 

339, 

.561057 

339. 

.580975 


407 67 6.0 










146; 8096 


1017; 4008 

; i 

H3SB 

9022 

1085; 4034 

m 

9069 


339.590961 I -8.636819 


-8.482400 


I BfcfclWAMW I 

KIlkEEl 


339.598123 | -8.477542 
66 

339.603519 I -8.465211 


-8.559631 


■BCEESBES 


339-625113 


339.627154 


339.628997 


339.641625 




339.651048 


339.655479 


339.656071 


339.661912 


■■ -m.-m .m 


339.676826 


339.681931 


339.682437 


339.684090 


339.688585 


339.697056 


339.702035 


339.704640 


339.704964 


339-705753 


ifcum.M.rea 


339.707062 


339.708279 


imiumn 


339.714394 


339.719062 


339.720638 


339.730426 


339.736522 




339.752527 


339.757206 


339.761885 


339.761725 


339.770798 


339.774642 


339.784544 


EiSSiSiiiii 


339.810757 



-8.455913 


-8.596113 


-8.645489 


-8.432215 
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Table 4. IDAS Asteroid Sightings 
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Table 5. Known Asteroids in the 15 June 1996 ISO Field 


(IDAS) Asteroid 

ID 


RSH 

M 

SNR 

UTC 








(D _ 

3 

339.57598 

-08.57781 

4032 

25.1 

06 58 

1999 AQ23 


339.57288 

-08.57761 




(1) 

1008 

339.58899 

-08.58097 

4339 

17.2 

09 25 

1999 AQ23 


339.58596 

-08.58042 











(2) 

4 

339.64707 

-08.46325 

5543 

22.0 

06 45 

17971 1999JZ50 


339.64713 

-08.46247 




(2) 

1003 

339.67256 

-08.45494 

5522 

24.3 

09 11 

17971 1999 JZ50 


339.67267 

-08.45481 





Table 6. Horizons Initial Heliocentric Osculating Elements for Ecliptic and 
Mean Equinox of J2000.0 and Epoch = 2001 -Oct-1 8.0000000 (TDB), for 1999 AQ23 
and 2001 -Apr-01 .0000000 (TDB), for 17971 1 999 JZ50 1 


Asteroid 

-IP . 

q 

e 

AP 


i 

H 

1999AQ23 

2451507.3543131 

2.304351999 

0.1217025 

97.89149 

135.42849 

14.72242 

13.5 

17971 1999JZ50 

2451472.0746345 

1.887653106 

0.168258407 

167.6773526 

128.5014528 

2.9364673 

14.8 


Tp = Time of perihelion passage (Julian date), q = Perihelion distance (AU), e = Eccentricity, AP = Argument of perihelion (deg), 
Omega = Longitude of the ascending node (deg), i = Inclination (deg), H = Absolute visual magnitude (G, the slope parameter, is 
assumed to be 0.15 for both asteroids). These data were obtained from Horizons on 9 Sep 2001 . 


Table 7. Aspect Data, Flux, and Derived Diameters and Albedos for Known 
Asteroids Potentially Associated With IDAS Sources 


Asteroid 

El 

HDist 

AU 

GDist 

AU 

EH 


D-STM 

km 

Ph-STM 

D-NEATM 

km 

Ph-NEATM 

1999AQ23 

18.2 

2.538 

2.053 

22.57 

4.08 + 0.15 

2.80 ±0.02 

0.90 ±0.01 

3.30 ± 0.02 

0.65 ± 0.01 

17971 1999 JZSO 

18.1 

1.892 

1.333 

31.01 

5.57 + 0.01 

1.50 ±0.01 

0.93 ±0.01 

1.73 ±0.01 

0.71 ±0.01 



Below 

are results assumin 

g the actual H is one full magnitude higher 1 

1999 AQ23 

19.2 

2.538 

2.053 

22.57 

4.08 ±0.15 

BsnnffiEi 

0.48 ±0.01 

2.95 ± 0.02 

0.32 ± 0.01 

17971 1999 JZ50 

19.1 

| 1.892 

1.333 

31.01 

5.57 ±0.01 

1.30 ±0.01 

0.50 ±0.01 

1.54 ±0.01 

0.36 ± 0.01 

Below are results assuming the actual H is 0.5 magnitude higher and the Flux a factor of 1 .5 greater. 1 

1999 AQ23 

18.7 

2.538 

2.053 

22.57 

6..12 ± 0.15 

3.00 ±0.02 

0.51 ± 0.01 

3.60 + 0.02 

0.34 ±0.01 

17971 1999 JZ50 

18.6 

1.892 

1.333 

31.01 

8.36 ± 0.01 

1.60 ±0.01 

0.52 ± 0.01 

1.90 ±0.01 

0.3710.01 
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Figure 1. Sample Region Coverage of a 17 x 17 Raster Map. Each V x V cell shows 
the total number of 20 second integration sets made on that area. 
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» v # 4 ♦ • ** « + .*** **** *# » * • * * % % *r • 

a 1 996 June 1 5 05:37:57 - 08:02:09 UT 



b 



1996 June 15 08:03:21 - 10:27:33 UT 



c 1997 June 10 00:19:15-02:13:35 UT 



d 1997 June 10 02:14:39 -04:08:59 UT 




e 



1997 June 10 09:14:03 - 11:38:14 UT 



f 1 997 June 1 0 1 1 :39:26 - 1 4:03:38 UT 


Figure 2. ISOCAM maps in ecliptic coordinates: 1996 June 15 (a, b) and 1997 June 10 (c, d e, f). North ecliptic latitude is up 
and east ecliptic longitude is to the right. Maps a, b, d, and f are 17x17 rasters while c and d are 15x15 rasters. 
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Figure 4. IDAS Inertial Point Sources. Squares outline the areas sampled nine times 
per map. The inertial sky within the union of the two large squares was sampled 36 
times (total integration time 720 sec per map pixel) and that within the two small 
squares 54 times (total integration time 1 ,080 sec per map pixel). The point size is 
proportional to the flux density, which ranges from 0.3 to 10.8 mJy. 
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Figure 5. IDAS Non-inertial Point Sources, i.e-., all sources from Fig 3 not plotted in Fig 
4. Sources from the 1 996 field are plotted in the top panel and those from the 1 997 field 
in the bottom panel. Squares outline the areas sampled nine times per map. The point 
size is proportional to the flux density, which ranges from 0.4 to 4.3 mJy. Circles 
indicate sources extracted from Maps 1 (top) or 3 (bottom), triangles those from Maps 2 
(top) or 4 (bottom), squares those from Map 5, and diamonds those from Map 6. 
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Figure 6. Observed (circles) versus predicted (triangles) positions for two bright known asteroids in the 
1996 ISO map. Each graph is 3’ on a side, the approximate size of the ISOCAM array. The small squares 
centered on the observed position are 6" on a side, the approximate size of the ISOCAM PFOV used. 







Correction Factor 


1.5 



Figure 7. Moving source flux correction. 
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These data provide the first direct measurement of the 12 pm sky-plane density for asteroids on the ecliptic 
equator. 

The median zodiacal foreground, as measured by ISOCAM during this survey, is found to be 22.1 ±1.5 
mJy per pixel, i.e., 26.2 ± 1 .7 MJy/sr. 

The results presented here imply that the actual number of kilometer-sized asteroids is significantly 
greater than previously believed and in reasonable agreement with the Statistical Asteroid Model (Tedesco, 
et al., 2002b, Astron J., to be submitted.). 
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